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Abstract 

Using a system composed of a Quadrupole and Ioffe Configuration (QUIC) trap loaded from a 
slowed atomic beam, we have performed experiments to observe the Bose-Einstein Condensation 
of Na atoms. In order to obtain the atomic distribution in the trap, we use an in situ out of 
resonance absorption image through a probe beam, to determine temperature and density. The 
phase space density (D) is calculated using the density profile and the temperature. We have 
followed D as a function of the final evaporation frequency. The results show that at 1.65 MHz we 
crossed the value for D expected to correspond to the critical point to start de Bose-Condensation 
of the sample. Due to the low number of atoms remaining in the trap at the critical point, the 
interaction produces minor effects and therefore an ideal gas model explains well the observations. 
We analyze the obtained low number in terms of efficiency of evaporation. The utility of an in situ 
detection is illustrated by measuring the harmonic gas pressure of the trapped gas in the route to 
condensation. 

PACS numbers: 03.75.Hh; 32.80.Pj 
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I. INTRODUCTION 



Bose-Einstein Condensation (BEC) is a field of research with vast interest mainly because 
it touches several subjects in physics. As a phase transition from a classical gas to a coherent 
matter-wave sample, it provides a large variety of interest in thermodynamic. As a heavy 
occupation of phase space, BEC is a nice object for quantum statistic. As a spontaneous 
symmetry breaking, it shows interest in quantum field theory. Many phenomena related 
to BEC in the past, like superfluidity and superconductivity can now be investigated in a 
controlled way. The field of atomic quantum optics has also presented many interesting 
possibilities related to BEC. Finally, BEC has opened up new windows into the quantum 
world in general. 

One of the interesting possibilities to be investigated using a BEC is the creation of 
condensate samples in a non ground state configuration. Such a non-equilibrium system can 
be formed if one first obtains the ground state condensate followed by a resonant pumping 
transfer to an excited state. Such schemes have been theoretically investigated by our group 

and shall be one of the goals of our experiments in the future. Another interesting 
possibility is the investigation of BEC from the thermodynamic point of view considering 
thermodynamic variables Q|. A sample of such investigation is presented in this paper. 

In this paper we report the first realization of BEC in Na atoms in our laboratory. We 
have followed, using absorption images, the temperature and density of atoms within the 
trap up to our optical resolution. Using these parameters we have followed the phase-space 
path that takes the sample to quantum degeneracy. We start presenting our experimental 
setup followed by the results and discussion. An application of our in situ observation is 
the determination of the thermodynamic parameters and evolution as the atoms are cooled 
to degeneracy, which is presented in the last section. 



II. EXPERIMENTAL SETUP 

In order to generate the conditions to attain BEC in Na, it is necessary to start with 
a thermal beam. We have loaded a Magnetic Optical Trap (MOT) from a slow atomic 
beam, followed by the necessary steps to BEC [3J. Our experimental system is presented in 
figure 1. An effusive sodium beam is decelerated in a tapered solenoid employing the dark- 



2 



FIG. 1: Schematic diagram of the experimental system. Atoms emerging from an effusive oven 
are decelerated and trapped in a MOT aligned with slowing tube. The magnetic field profile is 
presented in the figure detail, where z in given in centimeters and B in tesla. 

SPOT Zeeman tuning technique jj]. The process ends with a continuous flux of slow atoms 
emerging from the solenoid in the lower hyperfine state. Although the atoms are decelerated 
in the cycling transition 3S , i/ 2 (i r = 2) — > 3P 3 /2(F' = 3), as they approach the end of the 
slowing solenoid, the lower amplitude of the field as well as the configuration of the field 
lines create adequate conditions to optically pump the atoms to the 35*1/2 (F = 1) hyperfine 
ground state. At this point the deceleration process stops abruptly and the slowed atoms 
migrate out of the solenoid without interacting with the slowing laser, travelling about 40 
cm before being captured by the beams of a MOT. 

Three ring dye lasers (Coherent-699) provide the laser beams for the slowing, trapping 
and repumping transitions. All lasers are frequency stabilized and externally locked to the 
appropriated atomic transition, using vapor cell and saturated absorption signal. The slower 
laser is tuned to about 50 MHz to the red of the 3Si/%(F = 2) — > 3P 3 /2(F' = 3). The two 
other lasers are used to produce the trapping and repumping frequencies of the MOT. The 
trapping frequency is tuned -5 MHz from 3Si/2(F = 2) — > 3P 3 /2{F' = 3) transition while 
the repumping light is tuned to 3S , i/ 2 (F = 1) — > 3P 3 / 2 (i ?/ = 2) transition. 

The MOT operates as a dark-SPOT [5j, where the trapped atoms are confined in a 
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FIG. 2: Configuration of the QUIC coils of our experimental system. The coils for quadrupole field 
and the third coil for Ioffe is in conical shape to improve the optical access. 

hyperfine ground state which does not interact with the trapping laser frequencies. To create 
the dark- SPOT we placed a physical obstacle to block the central region of the repumper 
laser beam. In such scheme, the atoms located in the dark region are rapidly pumped to the 
lower energy ground state where they do not interact with the lasers, allowing to produce 
a dense sample of ground state atoms. This is a very important point to make an efficient 
transfer of atoms from the MOT phase to the magnetic trap. 

n 

The magnetic trap is the Quadrupole and Ioffe Configuration (QUIC) [6], whose coils ar- 
rangement is presented in figure 2. The system is composed of the two main coils (quadrupole 
coils) and an Ioffe coil. Each of the coils were winded using copper tube (with a diameter 
of 1/8 inch) and the number of turns in each coils is 50 for quadrupole and 30, for the Ioffe 
coil. For a quadrupole current of 230 A and the Ioffe of 220 A, the produced trap field in 
the Ioffe axis as well as the field contour lines in the plane x-z are shown in figure 3a and 
b, respectively. The QUIC trap confining potential is close to an harmonic oscillator, in our 
case, with the following frequencies: uj x = 2ir x 36.7 Hz; u y = 2tt x 120.8 Hz and uj z = 2n 
x 159.6 Hz. 

ri 

Once the atoms are in magnetic trap the evaporation process starts |7[. Through a one 
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FIG. 3: The produced trap field for a quadrupole current of 230 A and the Ioffe of 220 A. (a) Field 
profile in the Ioffe axis and (b) field contour lines on the plane x-z, in which each value correspond 
to the field in Gauss. 

photon RF transition 3Si/ 2 (F = l,m F = —1) — > 3Si/ 2 (F = l,m F = 0), the hottest 
atoms are transferred to the untrapped state carrying away energy. The RF determines the 
magnetic field at which position atoms undergo a spin-flip transition, been ejected from the 
trap. Thus the remaining atoms rethermalize through elastic collisions to lower tempera- 
tures. Therefore, sweeping the RF from a top value down, we force the evaporation until 
reaching the critical temperature to observe BEC. During the evaporation stage, the atoms 
must be shielded from any resonant light, that may cause optical transition destroying the 
sample. The RF antenna is a coil placed close to the cell having 5 turns of wire in a rect- 
angular frame of 2.5 x 7 cm. In our experiment we use two stages of evaporation: first on 
the quadrupole trap after compression and the final evaporation on harmonic profile of the 
QUIC when Ioffe coil has full current. 

Our transference procedure consists of loading the MOT with about 10 9 atoms and ini- 
tially transfer then to the quadrupole trap followed by compression during 2 seconds. After 
final compression in quadrupole trap, at a field gradient of 230 Gauss/cm, RF is applied to 
promote evaporation during 5 seconds from 100 to 30 MHz. Then after this first stage of 
evaporation, the Ioffe coil is turned on gradually during 2 seconds, producing the final field 
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FIG. 4: The full time sequence of the magnetic trap loading and imaging. 



configuration where the evaporation procedure is completed with RF from 30 MHz to final 
frequency (1.65 MHz to reach BEC). 

At different stages of the cooling process, either a time of flight (TOF) measurement can 
be performed to obtain the temperature or the in situ absorption imaging can be performed 
similarly to experiments already done [a]. Our imaging system, composed of lenses and a 
spatial filtered laser beam, allows us to obtain the absorption image with a resolution of 5 
fim. We have obtained absorption image with the magnetic field on. The image is normally 
obtained together with interferences and remaining defects of the optics. Those imperfections 
can be cleaned out using a Fourier Transform procedure to remove high frequency structures 
on the image. The full time sequence of the magnetic trap loading and evaporation as well 
as imaging timing is presented in figure 4. 



III. RESULTS AND DISCUSSION 



Before starting the presentation of our results we would like to present the overall sequence 
of absorption images given to the reader the general idea of the experiment. On figure 5 on 
the right hand side we show the typical size of the atomic cloud and the number of atoms. 
On the left hand side we have the typical temperature for each stage. In brief, the sample 
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FIG. 5: Sequence of absorption images given the general idea of our experiment. On the right 
hand side the typical size of the atomic cloud and the number of atoms. On the left hand side we 
have the typical temperature for each stage. 

starts with 10 9 atoms at 500 fiK and ends at 70 nK with a few thousand atoms in the 
condensate. 

When we start, the atoms are very much spreaded around the quadrupole trap. Starting 
the evaporation the Na atoms concentrate towards the center of the trap and the asymmetric 
geometry of the magnetic field shows up on the atomic cloud. Just after reaching BEC, the 
limit of resolution of our optical system appears to start to become important. 

These results were obtained using the relation between the images and the magnetic trap 
geometry. The QUIC potential can be expresses as 



Through absorption imaging one can obtain the values of a x and a z as well as the peak 
density at the trap center n , using a simple integration. Considering <r x and <j z one ob- 
tain the temperature of the cloud, since a hotter gas spread more in space according to 




(1) 



and the atomic density is found to be in good approximation with 
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FIG. 6: The peak phase space density as a function of final RF. For the final RF at 1.65 MHz, 
the Bose-Einstein Condensation is achieved. 



YT—mou^a^. Once we obtain no and T, the peak phase space density D is evaluated. Con- 
sidering D = noA 3 where A = ^ h mkT is the thermal de Broglie wavelength, as the peak 
phase space density, one obtain for Na that 

_ 4,84xlQ- 23 n 

u - t! W 

for no in cm -3 and T in Kelvin. 

According to previous work 9], once D reaches the value 2.612, the condensation starts 
to take place. For each final evaporation frequency we obtain D and the result is presented 
in figure 6 where the peak phase space density evolution can be followed and for the final 
RF frequency of 1.65 MHz, the Bose-Einstein Condensation is achieved. At this point we 
obtain D ~ 14, but still most part of the atoms are in the thermal cloud. Since at this 
point we have reached our resolution limit, thus we just overestimate the temperatures and 
underestimate the peak density, thus D must be larger than the reported value. We should 
observe that only a few thousand atoms are predicted to be in the condensate. 

According to characteristic values of our trap, the ground state size should have cr° ~ 
2.4 /mi, 0"° ~ 1.32 fim and a® ~ 1.14 /im, if interactions are not considered. Using those 
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values and the measured number of atoms in crossing the BEC line, we obtain the critical 
temperature Tc ~ 61 nK which is very close to the measured 70 nK. 

This kind of analysis is only valid if the system is close to an ideal gas. To verify if that 
is true we must compare the atomic energy of the ground state with the interaction energy 
at the critical point. The ground state is e$ ~ (wo = ^/uj^uj v uj 7 ), which correspond 

to 6.7 nK. On the other hand in the Thomas- Fermi approximation [lOj, the interaction 
energy/particle is given by Energy /particle ~ ^n JJ, where U = 4lT ^ a with a = 52(5) 
°o which is the scattering length. In the present case, Energy /particle ~2 nK and 
therefore we can neglect the interactions and consider the system as an ideal gas. 

The small number of atoms attained in our BEC can be justified through of an investi- 
gation of the variation of the number of atoms as a function of temperature. The obtained 
result is presented in the figure 7. As the evaporation on QUIC trap progress, two distinct 
regimes are observed. In a first regime NocT 137 , while for the second NocT a4 . Such de- 
pendences can explain the low achieved number. Efficient evaporation requires an elastic 
collision rate to be larger than the loss rate. For continued cooling, the so called runaway 
evaporation, the elastic collision rate must be maintained despite the atom number reduc- 
tion. The dependence of the atom number with temperature can be written as NocT s . If 
s>l, we have not reached the runaway evaporation condition and we lose more atoms than 
necessary. If s<l the condition for runaway evaporation is fulfilled, and better evaporative 
efficiency is obtained. 

In our case we observe that runaway evaporation is only obtained at the final stage of 
the evaporation process. Because of this, most part of the atoms are expelled from the trap 
remaining only the coldest part and this costs a large number of atoms, which explains the 
BEC with such small number. We have tried to optimize the first part of evaporation to 
improve the final result, however we have not yet been able to accomplish such improvement. 
Most probable reason is that for obtaining such improvement we have to slow down the RF 
sweeping speed, but in this case the lifetime of the trapped atoms due to background collision 
become a limiting fact on the process. Some modifications on our system must to be done 
to solve this problem. 
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FIG. 7: As the evaporation progress, two distinct regimes are observed. In a first regime NocT , 
while for the second NocT ' 4 

IV. EVALUATING THE HARMONIC PRESSURE DURING THE COOLING 
DOWN PROCESS 

Little attention has been paid to the thermodynamic of the cold trapped gas system. 
For this system, usual pressure and volume are no longer thermodynamic variables, since 
the confining potential is not homogeneous. Different from the rigid walls container, atoms 
trapped in a magnetic field interacts (with different amplitudes) everywhere in the potential. 
This results in a non homogeneous spatial distribution of gas. In this sense we can not 
carry on a thermodynamic analysis using conventional pressure and volume to investigate 
the confined gas thermodynamics. It is necessary to redefine the intensive and extensive 
variables. In a recent paper, V. Romero- Rochin has investigated a confined gas in a three 
dimensional harmonic potential and shown that pressure and volume are well replaced by 
new variables. The extensive variable "volume " V is now replaced by the inverse cube of the 
harmonic frequency (V = cUq 3 ). There is also a conjugated intensive variable, which plays the 
same role of usual pressure and is responsible for the mechanical equilibrium of the gas. Both 
variables work as "volume" and "pressure" in the thermodynamical sense. Considering 
a system composed of N particle in equilibrium at temperature T, the thermodynamic 
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treatment reveals that the state equation for such system, i. e., the combination of redefined 
V and P gives P): 

P= 2 -ujl(U ext ) (4) 

where (U ex t) represents an average integration over the external potential. This pressure is 
referred as harmonic pressure to distinguish the conventional pressure for gas in a box. In 
an explicit form, if the gas is subjected to the potential as Eq. (1) has a density profile 
n (x,y,z) as in Eq. (2) , one gets 

2 f 1 

P = 3^0 / dxdydz n (x,y,z) -m(uj 2 x x 2 + uj 2 y y 2 + u 2 z z 2 ) (5) 

In our experiment we have followed in situ n (x,y,z) and knowing well the potential 
parameters we have followed P as a function of T. Since J dxdydzn (x,y,z) =N one find that 
for the Gaussian density distribution, one obtain 



(u 2 y x +uy y +uy z ) (6) 



NC 

where <r x , a Y and a z depend on the gas temperature and C is equal to a dimensional numerical 
constant which depends on the single particle mass. 

In this work we have done the first set of measurements to prove Eq. (6) in a situation 
where u is kept constant. This is equivalent to an isocoric transformation of the trapped gas. 
As the evaporation takes place we record cx x and cx z and estimate <r y , by the relation between 
the oscillation frequencies. Thus with the total number of particle and gas temperature we 
create a plot of ^ as a function of T, as shown in figure 8, previously to reach BEC. 

The observed linearity of ^ with T shows that PV oc T, as expected since the limits 
presented by V. Romero- Ro chin Q in his work is the classical limit where the defined pressure 
and volume would obey the ideal gas state equation. We would like to explore the deviation 
of the state equation from the ideal gas case. That will allow us to finally evaluate the 
heat capacity and other important thermodynamic functions as the gas passes from the 
classical to the quantum regime. This work is now under preparation. For this purpose 
is extremely important the detection of the atomic distribution as the gas is still trapped 
(in situ detection). A full investigation of harmonic pressure under isothermal or isochoric 
transformation will be implemented in our laboratory. 
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FIG. 8: Graphic showing the evolution of normalized harmonic pressure as a function of 
temperature T. 

V. CONCLUSIONS 

In this work we demonstrate that we have achieved Bose-Einstein condensation in a 
sample of Na atoms by in situ observation of the atomic density profile and temperature, 
resulting in the evaluation of the phase-space density. The overall analysis allows us to 
understand that the final low number of atoms was consequence on the fact that runaway 
evaporation only took place at the final stage of the process. The route to BEC was also 
analyzed on a thermodynamic view showing the evolution of the gas pressure as the temper- 
ature was lowered. The in situ observation of the gas will allow us further thermodynamical 
analysis of the system as it approaches BEC. 
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